Experimental signatures of cosmological neutrino condensation 
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Superfluid condensation of neutrinos of cosmological origin at a low enough temperature can 
provide simple and elegant solution to the problems of neutrino oscillations and the accelerated 
expansion of the universe. It would give rise to a late time cosmological constant of small magnitude 
and also generate tiny Majorana masses for the neutrinos as observed from their flavor oscillations. 
We show that carefully prepared beta decay experiments in the laboratory would carry signatures 
of such a condensation, and thus, it would be possible to either establish or rule out neutrino 
condensation of cosmological scale in laboratory experiments. 
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The flavour oscillations of neutrinos and the acceler- 
ated expansion of the universe are two fundamental prob- 
lems in particle physics and cosmology. The neutrino os- 
cillations indicate that neutrinos are massive and their 
masses are slightly different for the different flavors. On 
the other hand the accelerated expansion of the universe 
requires an existence of some form of dark energy or a 
cosmological constant. Due to the proximity of the dark- 
energy scale with the neutrino masses, there have been 
a lot of efforts in finding an unified solution to both the 
problems 

Recently, there has been a lot of interests in neutrino 
condensation on the cosmological scales Q. It has been 
argued that the background neutrinos under certain con- 
ditions can become superfluid [3, |j] . The neutrino con- 
densates can then contribute towards the accelerated ex- 
pansion of the universe and they can also generate 
the Majorana mass term. In this work we investigate this 
scenario further with the aim of investigating experimen- 
tal signatures of the background neutrino-condensates. 

Our arguments are based on the phenomenon called 
the "Andreev reflection" which is extensively studied in 
the literatures of superconductivity 0] ■ Before explaining 
the Andreev reflection, we first mention the general ideas 
for the neutrino condensation. Our analysis is based on a 
very simple and elegant model given by Caldi and Cho- 
dos [3j, although the results are more general and can 
be applied to the other mechanisms also. This picture 
requires that the background neutrinos to be in a de- 
generate gas/liquid state. The idea of degenerate neu- 
trinos in the universe is not inconsistent with the stan- 
dard model of cosmology It is well known that at a 
low enough temperature a gas of degenerate fermions un- 
dergo a phase-transition to the superfluid state, if there 
is an attractive-interaction between the fermions in any 
angular- momentum channel [10l| . However such transition 
for the Dirac neutrinos within the standard electroweak 
interaction may not be possible (3 (see also Ref.[ll|). In- 
terestingly, they also suggest that the degenerate Dirac- 
neutrinos in the presence of some new (unknown) attrac- 
tive interaction, howsoever weak, can form condensates. 

The condensation can be studied by using the La- 



grangian density for the four-fermion interaction within 
the the mean- field theory frame work 0, 0, [ll|. From 

this one can find the BCS-type superfluid condensation 

i 

with the gap parameter given by A ~ fj,e v 2 g 2 f where 
\x is the neutrino chemical potential and G signifies the 
strength of the four-fermions coupling. Such condensates 
can naturally give rise to Majorana mass of the neutrinos. 
In fact for the values of fi allowed by the big-bang nucle- 
osynthesis, it is always possible to choose the values of G 
such that the cosmological constant, A ~ G 2 | < vv > | 2 , 
is of the order (10 _3 eV) 4 = (meV) 4 and the value of neu- 
trino mass, m v ~ A = G | < vv > |, can be anywhere 
between a meV and few hundred MeV . Neutrino mass 
generated through such a mechanism will not have any 
significant contribution to the dark matter in the uni- 
verse. However, if we allow the chemical potentials for 
different neutrino species to vary, the condensates could 
depend non-trivially on flavor, leading to an interesting 
spectrum of neutrino masses and mixing. This suggests 
a possible mechanism for neutrino oscillations without 
violation of the lepton number conservation at the mi- 
croscopic level. 



ANDREEV REFLECTION 

It is not difficult to understand the Andreev reflection 
phenomenon qualitatively. Generally the quasi-particle 
states propagating in a superfluid with energy e below 
the gap energy A are forbidden. But the sub-A transfer 
is possible if an incoming electron with e is transfered 
along with an another electron by forming a Cooper pair 
into the superconductor. In terms of single particle state 
this can be described by the reflection of a hole in the 
normal metal. Thus the Andreev reflection of a hole (or 
an electron) is equivalent to the transfer of a Cooper 
pair in (or out of) the superconductor. But the single 
electron transfer with the sub-gap energy is not allowed 
in the superconductor. Quasi-particles of energy e > A, 
propagate in the condensate as a massive particles, of 
mass equal to A [7|. 

In the mean field description, the starting point is an 
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effective Hamiltonian with four-fermion interaction given 
by, 

H= f d 3 x[j2 E *)Mr,t) 
o=t,; 

+A(r, *)V>K r ' *) + A*(r, *)^(r, ij^ffc *)] W 

where, Ep is the Fermi energy (~ chemical potential), 
and the "order parameter" field A is 

A Q(3 (r, t) =< ijj a {r, t)ip p (r, t) >= A(r, i)e a/9 (2) 

The equation of motion obtained from the Hamiltonian 
above is known as Bogoliubov- de Gennes equation in the 
theory of superconductivity/ superfluidity and is given 
by, 
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(3) 



These linear equations can be treated as equations for 
single particle wave functions instead of operators. For 
A = 0, the equations decouple and we have the plane 
wave solutions, with the dispersion relation, E = e q 
for particles and E = —e q for holes (antiparticles), 
e q = q 2 /2m — Ep. This is the case with the normal con- 
ductors (metals). When A = const. ^ 0, we still have 
plane waves as solutions, however, with the change that 
quasiparticle spectrum within the condensate is gapped, 



E 2 = e 2 q 



(4) 



The quasiparticle wave functions are, 
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exp(tq- ■ r — lEt) (5) 



junction of the materials is located at z = 0. At the con- 
ductor side (z < 0), A = while at the superconductor 
side (z > 0), A = constant. Let us now suppose that 
particles of energy, E, from conductor side of the slab, 
hits the junction. To see the details of what happens 
at the boundary, we need to solve the stationary Bogoli- 
ubov - de Gennes equation, Eq.(3), with the boundary 
conditions : (1) for z — > — oo excitations are particles or 
holes (conductors), (2) for z —> oo excitations are quasi- 
particles (superconductor), (3) for z = wave functions 
and their first spatial derivatives are continuous at the 
boundary at z = 0. It is not hard to solve this prob- 
lem and detailed solutions can be found in reference Q- 
The net result is that the probability current through the 
junction is, 







for E < A 
for E > A 



(6) 



where vp — y / '2Ep/ni, is the Fermi velocity. The reflec- 
tion and transmission coefficients are, 



Rhole — 1. 
u — E-l 



T„ 



T 
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for E < A 
for E > A 



(7) 



These equations clearly show that an incoming particle 
hitting the junction from the conductor side with energy 
below the gap can not cross the junction, instead a a hole 
is reflected back in the opposite direction. This is not 
difficult to understand physically- near the boundary, the 
incoming particle creates a particle and a hole of equal 
and opposite momenta from the Fermi sea, pairs up 
with the particle and disappears in the condensate and 
the hole is reflected back. This is the Andreev reflection 
phenomenon. For E » A, R hoie -> 0,T guos , -> 1. For 
energy, E, not too high above the gap, A, there is partial 
transmission of quasiparticles and partial reflection of 
holes. 



where A ± = ^{l±£/E), £ = - |A| 2 , 4/2m = 
Ep ± £, 8 is the phase of the complex order parameter, 
and D and F are some constants. The wave functions, 
<j) p and 4>h correspond to particle and hole like excita- 
tions of mass A. However, note that they are not in 
the eigenstate of particle number or electric charge. For 
E >> |A|, they reduce to wave functions for particles 
and holes respectively. When the quasi-particle excita- 
tion energy, E < |A|, the momenta are complex, and 
therefore, these modes do not propagate in the conden- 
sate medium. When such excitations reach the normal 
conductor-condensate junction from the conductor side 
of the arrangement, we have, what is known as Andreev 
reflection in the theory of superconductivity. 
Let us consider a normal conductor-superconductor slab 
aligned along the z — axis. Assume r = (0,0, z). The 



EXPERIMENTAL SIGNATURES 

To understand the role of Andreev reflection in the 
context of neutrinos, we have to find the fate of the con- 
densate within materials. Cosmological neutrinos are of 
low energy and large wavelength, and threfore, the size of 
the cooper pairs in the condensate would be large. Weak 
interaction of neutrinos with matter is mediated through 
charge currents as well as neutral currents. However, the 
cross section of such interactions is extremely small. For 



E v = 1 M eV, this cross section is ~ 10 



and for 



E u = 1 eV, it is 
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-56, 



Thus it would seem that 



matter is almost transparent for neutrinos. However, it 
turns out that for low energy neutrinos there is coherent 
neutrino processes in matter 12| . For momentum ~ a 
few hundred MeV/c, neutrinos would elastically scatter 
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from nuclei with cross section proportional to A 2 . For 
lower momenta it will coherently scatter off the entire 
atom, meaning all the nucleons and electrons in the atom. 
Subsequent lowering of momenta involves, coherent scat- 
tering of neutrinos from all the atoms contained in the 
wavelength. Thus the scattering cross section would grow 
as the square of number elementary particles ( quarks & 
lcptons) contained within the wavelength. Thus, it would 
seem that the cooper pairs could be broken and conden- 
sate could be disordered by coherent scattering processes. 
We know from work of Caldi and Chodos [3| that the 
value of the order parameter is decided by the chemical 
potential and the coupling constant of the assumed " new 
interaction" . From the big bang nucleosynthesis argu- 
ments, we get only the upper bound on the combined 
value of the chemical potentials of all the three flavors 
of neutrinos. We may get an estimate of the range of 
the "new coupling" from the formula for the neutrino 
mass and the cosmological constant obtained from con- 
densate order parameter. The neutrino mass is given by, 
m v ~ A = G' 2 \ < vv > | and cosmological constant by, 
A — G 2 \ < vv > | 2 , and therefore, G 2 ~ m 2 / A. If 
m„ ~ meV (~ eV) and A ~ (meV) 4 , range of the 
coupling G, is of the order of meV. But for vn v ~ eV, 
the range changes to micro electron volts. Moreover, for 
different flavor of neutrinos, we expect that the chemical 
potentials as well as the coupling constants are slightly 
different. Thus there is large uncertinity in the values of 
the relevant parameters and it is hard to decide whether 
the cooper pairs are stable or disordered by coherent scat- 
tering processes within ordinary materials. Therefore, for 
experimental signatures of condensates, we will consider 
both the possibilities. In the first case, the cooper pairs 
are assumed to be broken within ordinary materials and 
the material surface forms boundary with the cosmolog- 
ical condensate. In the second case, we assume that the 
condensate penetrates ordinary matter and the cooper 
pairs are intact. 

We have, so far, considered only the neutrinos of cosmo- 
logical origin. However, there exist many other sources of 
neutrinos in nature as well laboratory. The energy spec- 
trum of neutrinos produced in these processes is very 
wide, but most of the flux is concentrated in energy of 
the order of a few MeV. With the cosmological context in 
mind, we are interested in properties of low energy neu- 
trinos. It turns out that the radioactive decay of tritium 
contains low energy flux of antineutrinos. We will, there- 
fore, consider an experimental set up with radioactive 
tritium sample undergoing beta decay. The experimental 
set up should be conceptually similar to the arrangement 
used in the KATRIN experiment The tritium nu- 

cleus decays to, 3 Hi — > 3 He2 + e~ + v e , with half life 
of Ti/2 ~ 12.3 years. The total energy released in the 
decay is E = 18.6 KeV. Beta deca y e nergy spectrum 
is analyzed by using the Kurie plot [141 ] ■ The electron 



spectrum in the allowed beta decay is, 

N e {E e )dE e cx F(Z,E e )y/E 2 ~ m 2 e E e (E a - E e ) 

x V(£o - E e ) 2 - ml dE e (8) 

where F(Z,E e ) is the known Coulomb factor, Eq is the 
total energy released in the beta decay of the nucleus, 
E e , is the energy carried away by the electrons. The 
plot of the Kurie function, [N e {E e )/{F{Z, E e )p e E e )] 1 / 2 
versus electron energy E e should be a straight line 
when rriy — but should be of a different shape near 
the end point when m v ^ 0. We are interested in 
the properties of the neutrinos released in the beta 
deacy assuming the presence of cosmological background 
neutrino condensation. For clarity and simplicity, we 
assume that the background cosmological condensate is 
made from cooper pairs of antineutrinos. For neutrino 
condensates, one can use essentially the same arguments 
in a slightly different manner. In this situation, we have 
the following two possiblities: 

(1). The cooper pairs of the cosmological condensate 
are formed by antineutrinos of low energy and large 
wavelength. The pair binding is not too strong even 
by the standard of weak interactions, and therefore, 
within the radioactive material (tritium) the condensate 
is disordered by the process of coherent scattering. The 
condensate outside the material remains unchanged. 
We have antineutrinos both from the disordered cos- 
mological condensate and the the radioactive beta 
decay within the tritium sample. We encounter here 
a situation similar to one in condensed matter system 
with metal-condensate boundary. We, therefore, refer 
to Andreev reflection phenomenon to understand the 
physical processes in the material, condensate and their 
boundary. Antineutrinos ( and neutrinos if any) of 
energy lower than the value of the condensate order 
parameter |A|, are trapped inside the material. If any of 
these sub-barrier antineutrinos, of momentum p, reaches 
the condensate boundary, a neutrino-antineutrino pair of 
equal and opposite momentum, — p and p, is created at 
the junction. The newly created antineutrino combines 
with the old incoming antineutrino, forms a cooper 
pair and disappears in the condensate. The neutrino 
(hole) is reflected back. Similar phenomenon takes 
place when the reflected neutrino reaches the boundary. 
However, those with energy greater than |A|, move into 
the condensate and propagate as quasiparticles of mass 
|A| , their dispersion relation being, E = ^/p 2 + |A| 2 . 
Thus the energy of the antineutrinos from the beta 
decay in the sample covers the entire spectrum of the 
released energy and thus the Kurie plot in this case is a 
straight line (as in the case of massless neutrinos) with 
end point at the maximum energy, E — 18.6KeV (with 
a few extra events near the end point which will be 
explained latter). Such a mechanism of mass generation 
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does not lead to violation of lepton number at the 
microscopic level, and therefore, also implies that there 
is no neutrinoless double beta decay flBl ] . 
Andreev reflection process creates a mixture of sub- 
barrier energy neutrinos (holes ) and antineutrinos 
within the sample. The process 3 -Hi + v — > ^He^ + e~ is 
energetically feasible and has the same cross-section as 
3 Hi — > 3 He2 + e~ + v. Reaction cross section for such 
a processes involving low energy neutrinos is very small. 
However, in an experiment run over very long time some 
low energy antineutrino excess event, which would show 
up as excess electrons near the end point of Kurie plot, 
should be seen. The estimated number of such electrons 
is approximately 50% more than what is expected 
from the standard model. From the end point of the 
spectrum where the excess events start appearing, one 
can estimate the value of the order parameter. It should 
be noted that the background cosmological neutrinos, 
in the absence of condensation would also create excess 
events near the end of of the beta deacy spectrum, 
however unlike in our case, the energy spectrum will 
have a monoenergetic peak at E — Eq + m v , where Eq is 
end point of spectrum and m v is the neutrino mass [16| . 
Experiments such as KATRIN is designed to search for 
beta decay to low energy antineutrinos and it is in these 
experiments that one would expect such events to be 
established or ruled out. 

(2). The cosmological condensate penetrates the tritium 
sample and there is no change either in the supcrfluid 
ordering or the value of the order parameter. As ex- 
plained earlier, single particle states of energy less than 
the value of the order parameter can not be sustained in 
the condensate. Because of this constraint on the avail- 
able phase space, radioactive tritium can not decay into 
sub-barrier energy antineutrinos. However, it can decay 
into antineutrinos of energy greater than the value of the 
order parameter, |A|. These supra-barrier neutrinos pro- 
duced in the beta decay in the nucleus have wavelength 
shorter than the neutrinos in the cooper pair and they are 
in eigenstate state of particle number operator within the 
cooper pair length scale. Therefore, in the calculation of 
beta decay amplitude they are considered massless, with 
dispersion relation, E v — cp v (elsewhere, c = 1). In the 
derivation of Kurie function in this case, the important 
point that we have to keep in mind is that E u > A, 
because there is no single particle state below this en- 
ergy. Thus the Kurie plot is still a straight-line but the 
end point has a cut off at energy equal to Eq — |A|. How- 
ever, within the condensate these antineutrinos will prop- 
agate as quasi-particles (in Bogoluobov state) of mass 
|A|, with dispersion relation given by, E — y/p 2 + |A| 2 . 
As in the previous case, the the antineutrinos (and neu- 
trinos) propagate as massive quasi-particles having Majo- 
rana mass but there is no violation of the lepton number 
conservation at the microscopic level, and therefore, there 



is no possibility of neutrinoless double beta decay. The 
massive quasi-particles are not in the eigenstate states of 
particle number or the electronic lepton number opera- 
tors. This is also the necessary requirement for neutrino 
oscillations. Thus, slightly deferent values of the conden- 
sate order parameters for three flavors of neutrinos would 
nicely account for neutrino oscillations. 
In both the cases discussed above, we find that the value 
of the order parameter (and thus the mass of the neutri- 
nos) can be estimated from the end point behaviour of 
the beta decay spectrum. The rate of neutrinoless dou- 
ble beta decay is known to be proportional to the square 
of neutrino mass. Thus, we can estimate this decay rate 
and expect that in a suitably designed experiment, it can 
either be established or ruled out. 



CONCLUSIONS 

We have discussed in some details how the superfluid 
condensation of background neutrinos of cosmological 
origin can generate neutrino mass for flavor oscillations 
and also account for the accelarated expansion of the uni- 
verse. Such a mechanism does not seem to be in conflict 
with any aspect of the standard model of particle physics. 
However, it indicates the existence of new physics in the 
form of new attractive interaction among the neutrinos. 
We have shown that, in carefully prepared beta decay 
experiments, this model would lead to experimental sig- 
natures as listed below: 

(i) . Absence of neutrinoless double beta decay is a nec- 
essary requirement. If neutrinoless double beta decay is 
observed, the mechanism of mass generation and oscilla- 
tion discussed in this paper should be considered redun- 
dant. 

(ii) . The Kurie plot is a straight line with either of the 
two possibilities, (a) the end point of the plot extends up 
to the total beta decay energy released, Eq = 18.6KeV , 
with a few extra low energy events in an experiment run 
over sufficiently long time (such events are known to ex- 
ist, however, they are still within the experimental error 
bars and can be settled only by future experiments fl^j). 
(b) the straight line terminates at Eq — |A|, where |A| is 
the superfluid order parameter. 

Thus, it should be possible to either establish or rule out 
cosmological neutrino condensation from carefully pre- 
pared beta decay experiments in laboratory. 
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